A Droplet within the Spherical Model. 
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Abstract. Various substances in the liquid state tend to form droplets. In this 
paper the shape of such droplets is investigated within the spherical model of a 
lattice gas. We show that in this case the droplet boundary is always diffusive, 
as opposed to sharp, and find the corresponding density profiles (droplet shapes). 
Translation-invariant versions of the spherical model do not fix the spatial lo- 
cation of the droplet, hence lead to mixed phases. To obtain pure macroscopic 
states (which describe localized droplets) we use generalized quasi-averaging. Con- 
ventional quasi-averaging deforms droplets and, hence, can not be used for this 
purpose. On the contrary, application of the generalized method of quasi-averages 
yields droplet shapes which do not depend on the magnitude of the applied exter- 
nal field. 
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1 Introduction. 

The purpose of the Gibbs distribution and the statistical physics in general is to 
provide a bridge between microscopic and macroscopic phenomena. It would be a 
mistake though to think that the bridge is in any sense similar to one of, for instance, 
London bridges. The difficulties of actual "traveling" through that bridge in the case 
of finite-dimensional models (with d > 2) are so formidable, that a better name for 
the bridge might have been a labyrinth. Therefore it was not really surprising that it 
were mathematicians (not physicists) who actually managed to walk through a few 
bridges connecting microscopic interparticle interactions with macroscopic shapes of 
droplets of condensed matter. 

In a colossal effort beginning from works of Minlos and Sinai [8j |9] it was shown 
that for sufficiently low temperatures typical configurations of discrete lattice models 
look like a macroscopic droplet of one phase surrounded by another phase. The 
droplet border is sharp in the macroscopic scale, that is, the magnitude of typical 
fluctuations of the corresponding long contour is much smaller than the linear size 
of the droplet. Somewhat later the exact shape of the droplet in the case of the 
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Ising model was also found. It is given by a curve minimizing the corresponding 
Wulff functional. For an accurate account of all "twists and turns", which should 
be dealt with in order to walk through the bridge-labyrinth connecting the shape of 
a macroscopic droplet with the interaction of Ising particles, see the books [U [12] • 

The behaviour of continuous lattice systems is different from that of the above 
mentioned discrete systems, and properties of Gibbs states of, for instance, O (re- 
models are much less understood than those of various discrete models. Fortunately, 
there is a continuous lattice model which macroscopic properties can be derived from 
the interaction of microscopic variables with only modest efforts. This model is the, 
so-called, spherical lattice gas [6j [11] . 

The authors of the paper [11] studied the equation of state of the model — the 
behaviour of pressure as a function of specific volume and temperature. At the time, 
derivation of equation of state from a microscopic interaction was an achievement 
on its own right. Therefore, it is not surprising that such a question as the shape 
of a droplet of condensed "spherical matter" was not even considered. Moreover, 
the spherical lattice gas with cyclic boundary conditions is a translation invariant 
model. Hence, the center of the droplet is uniformly distributed over the available 
volume, and the constant average values of microscopic variables do not reveal the 
droplet shape. 

At nearly the same time, in the paper [1] similar phenomena related to invari- 
ance of correlation functions were termed the degeneracy of equilibrium state, and 
a heuristic procedure to remedy the situation — the quasi-averages — was adver- 
tised. Calculation of quasi-averages involves switching on an appropriate symmetry- 
breaking field of a magnitude proportional to e, and switching off the field by sending 
e 4- after the thermodynamic limit. It is shown in the present paper that the con- 
ventional quasi-averages allow one to find the true shape of the droplet only in 
a fortuitous situation when the symmetry-breaking field is already similar to the 
droplet shape we are looking for. The field of a different shape not only fixes the 
location of the droplet but also deforms it unrecognizably. Nevertheless, one can say 
that certain tools for finding the shape of a droplet in the spherical lattice gas were 
available already in early sixties, but investigation of that kind can not be found in 
the literature. 

Two decades later the dominant terminology became pure phases and mixed 
phases (instead of non-degenerate and degenerate equilibrium states, respectively), 
see, e.g., the book [7]. More importantly, a simple criterion allowing one to answer 
the question whether a phase is pure or mixed became widely known: a phase is 
pure if the covariance of microscopic dynamic variables associated with nodes j 
and k tend to zero as the distance \j — k\ increases. Pure phases can be obtained 
with a help of quasi-averages, or by using appropriate boundary conditions, or by 
calculating conditional distributions. Mixed phases can always be represented as 
linear combinations of pure phases. 

Gersch and Berlin, see [6], found the (constant) expected values, covariances, and 
conditional expected values of microscopic variables of the spherical lattice gas. The 
derived covariances do not tend to with the distance between the corresponding 
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microscopic degrees of freedom, hence, the natural phase of the lattice gas is not 
pure. In order to obtain the droplet shape in a situation like that, it is necessary 
to decompose the mixed phase into pure components. However, apparently, the 
authors of the paper [5] did not realize that. 

In the present paper we extract pure phases of the spherical lattice gas using 
(generalized) quasi-averages, see [21 EJ- The investigation of the properties of the 
pure phases shows that the droplet in the spherical lattice gas is always diffusive. 
That is, the boundaries of the droplet are not sharp, not even in the macroscopic 
scale. The constant levels of the expected values of microscopic variables look like 
rounded squares, although not exactly the same ones as the rounded squares de- 
scribing the sharp boundaries of the droplet within the 2D Ising model of a lattice 
gas. 

The rest of the paper is organized as follows. Section 2 contains a precise defini- 
tion of the spherical lattice gas. It also contains some well known technical results 
for the use in the later sections. Section 3 introduces macro states and summarizes 
the main results of the paper. In Section 4 we calculate the distributions of micro- 
scopic random variables in the mixed phase of the spherical lattice gas. In Section 5 
we use the Lagrange method to find pure macroscopic phases for the lattice gas with 
periodic boundary conditions at zero temperature. Section 6 is the main part of this 
paper. There we use the method of quasi-averages for extracting pure macroscopic 
phases. The results of the paper are discussed in Section 7. 

2 The model and useful facts. 

The spherical lattice g cLS IS cL collection of random variables {xj,j G Z d } placed at 
sites of an integer (^-dimensional lattice, Z d . Every site j G Z d is specified by its d 
integer coordinates (ji, J2, ■■■,jd)- I n the present paper we consider the case d > 3. 

To define the distribution of random variables at all sites of the lattice, we first 
specify the joint distribution for the random variables in a finite cube 

V n = {j eZ d :l < 3u <n,u = 1,2,..., d} 

containing N = n d sites, and then pass to the limit n — » oo. To avoid unnecessary 
complications we impose periodic boundary conditions in all dimensions. 

It is instructive to consider also a stretched model defined on the parallelepipeds 

T n = {j G Z d : 1 < j v <n v ,v = 1, 2, d}, (1) 

where n\ — (1 + 5)n, 5 > 0, and n u = n, for u = 2,3, d. Save for one side being 
longer than the others, the definition of the stretched model is exactly the same as 
that of the conventional spherical lattice gas. 

The Hamiltonian. 
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The random variables located in the rectangle V n interact with each other and with 
the external field {hj,j G Z d } via the Hamiltonian 



H n J y ' TjkXjXk y ' hjXj, 
jMVn jev n 



(2) 



where J > 0, and Tjk are the elements of the nearest-neighbour interaction matrix. 
In this paper the field {hj,j G Z d } is used as a technical tool, so that, it should 
not necessarily be physically sensible. We let its magnitude to depend on the size 
of rectangle V n : 

h =n^b„ jeZ d , (3) 
where the absolute values of bj are bounded by an independent of n constant b. 

The interaction matrix. 

The elements of the interaction matrix T describe the usual nearest neighbour in- 
teraction on a square lattice, and they are given by 



ze{i,2,...,4V 



where 

«*■»-{£ V&X, 

is the Kronecker delta. 

The coefficients J^(j u , k u ), for v = 1,2, ... ,d, are the elements of the matrix 
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The eigenvalues and orthonormal eigenvectors of the matrix J are given by 

2n(l - 1) 

A; = cos , / = 1,2, ... ,n, 



n 



and 



u w = < -u^/ = \ / — cos 



2n(l- l)(m- 1) 7T 
n 4 



, Z = l,2,...,ra. 



m=l 
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The eigenvalues of the interaction matrix T are the sums of eigenvalues of the 
matrix J: 

d 

At = ^\ > k= (k u k 2 ,... : k d ) eV n . 

The corresponding orthonormal eigenvectors are the products of eigenvectors of the 
matrix J: 

™ {k) = 1^ = {\uf; ] \ , k = {k u k 2 ,...,k d )eV n . (4) 

Note that the second-largest eigenvalue A s = d— 1 + cos(27r/ri) of the interaction 
matrix (which will play an important role below) is 2d times degenerate. At the 
same time, in the case of the stretched model the second-largest eigenvalue A s = 



d — 1 + cos 



2tt 



is only twice degenerate. 



_(l+8)n 

The Gibbs distribution. 

The joint distribution of the random variables {xj,j G V n } is specified by the 
usual Gibbs density 



p({xj,j G V n }) 



e -IBH n 



with respect to the "a priori" measure 

pt n (dx) = 5 I J2 x j - P N I 6 I E x ) - N ) II dx r 
The first delta function imposes the typical for gas models density constraint 

E ^ = p. 

the second one imposes the usual spherical constraint 

jev n 

The normalization factor (partition function) Q n (p) is given by 

/OO fQO 
... e-P H "fi n (dx). (5) 
-oo J —oo 



3 Macro states and the main results. 

The usual Gibbs states provide a detailed microscopic description of a thermody- 
namic system. In some sense the amount of available detail is too big: there can be 
a macroscopically inhomogeneous structure in the thermodynamic system, but its 
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shape can not be described within the Gibbs-state framework. In this situation an 
introduction of a rougher (reduced) description seems justified and useful. 

We define macro states as the following continuum limit of the original lattice 
system. The limiting configurations are realizations of random functions defined on 
the (i-dimensional rectangle [0, l] d : 



For any 7 e [0, l] d the random variable £(7) is defined as the following limit in 
distribution 



where [y], is the integer part of y. Of course, for correctness of the above definition 
we need some kind of continuity in the system, so that for two sequences j(n) and 
k(n) with the same limits 

Jhn (ji W/™, • • ',jd(n)/n) = Km (A* (n)/n, . . . , k d (n)/n) = (71,72, ... ,7«0> 

we also have identical limits of the corresponding sequences of random variables 
{a^(n)K£Li an d {%k(n)}%Li- Although in discrete spin systems that kind of continuity 
may well be missing, we do not worry too much about that, because it should be 
possible to surpass this technical problem one way or another. 

Thermodynamic random variables £(7) and x(S) are limits of the random se- 
quences a% 1 n],[72n],...,[ 7d n]) and x (lSin},[5 2 n},...,[5 d n}) separated by a distance of order n. 
Hence, in the continuum limit the random variables a; (7) and x(8) with 7 7^ 5 are 
independent due to the exponential/power-law decay of correlations in pure phases 
of high/low temperature regions. Therefore, a pure macro state is completely char- 
acterized by individual distributions 



Macro states defined above contain complete information about the microscopic 
individual distributions of the random variables Xj. It might be desirable to achieve 
further reduction of description by using, for instance, Kadanoff 's blocks in the defi- 
nition of the continuum limit. That should reduce the range of limiting distribution 
outside the critical lines/points to just the normal distribution. Even on critical lines 
the asymptotic distributions of large Kadanoff's blocks should be of quite limited 
variety. 

The main results of the present paper obtained for the low-temperature region 
j3 > (3 c {p), see Eq. (jHJ), can be stated as follows. 

1. The expected values of the random variables Xj in the natural state of the 
spherical lattice gas are translation invariant and equal to the density, (xj) = p. 
Their variances and covariances in the limit n — > 00 are given by 



{^(7)} 7 G[0,l]d = { x (li, 72, • • • , Id)} 



71.T2,— >7d6[0,l]- 



^^ X {hin],b2n],...,[y d n})i 



= iSL Pr « X (hxn],h2n],...,[ldn]) < V ■ 



Var(xj) 



W d (d) 1-p 2 



2/3J 2 



( 
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Cav(x j,xi) = c(j, + (l- P 2 ) 1- 



/3 c \lA 2n(j v -l v ) 

— - > cos 

P J df^l n 



oil) 



where c(j, Z) are the covariances of the microscopic variables in the ordinary 
spherical model, see Eq. f l2Tj) . 

The random variables Xj within the stretched model admit the following rep- 
resentation 



\ 



2(1 -P 2 



sin(27rC7)+^ U 



W d {d) ~ 
2(5J 



where the random variable U is uniformly distributed on the interval [0,1], 
J\fj(a, b) are normal random variables with mean a and variance b, and 

Cov(Afj,M) = c(j, I) -> as dist(j, Z) -> oo. 

3. For T = the pure equilibrium states of the model (ground states) are given 
by 

2tt(> - 1) 



a. 



P+v /2(l-p2)^ ri/COS 

v t^i V n J 

where a u and r u are any constants satisfying a u £ [0,2tc] and I]^ =1 r 2 = I. 

In the presence of a symmetry-breaking field 

if h 



en 
0, 



1. 



otherwise, 

the expected values of the random variables Xj are given by 



( x o) 



en 



(1-0/3 



-.71 



/ 2C , n -(l+M)/3 

2Jv^C 

if < 5 < 1 and ji < n/2. In this case all condensed "spherical matter" 
gathers in a narrow strip around the plane where the field is applied. The 
shape of the droplet is clearly deformed by the field. 

If 5 = 1 the distributions of the random variables Xj still depend on the 
magnitude of the symmetry breaking field, but in the limit n — > oo they remain 
proper. In particular, the finite expected values of the random variables Xj are 
given by Eq. (|2T)|) . The droplet shape is still deformed by the field. 

If 1 < 5 < d— 1, then the limiting distributions of the random variables Xj do 
not depend on the magnitude of the symmetry breaking field. In particular, 
the expected values of the random variables Xj are given by 



( Xj ) = p + 



\ 



2(1 



1 



P 



COs(27T7i) + o(l). 



The droplet shape is determined by the field (for different field types we obtain 
different shapes), but it is no longer deformed by the field (does not depend 
on varepsilon). 
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4 The mixed phase. 



Main thermodynamic properties of the spherical lattice gas with cyclic boundary 
conditions were derived in the papers [H [11] . Although the free energy density of 
the model is not sensitive to the type of boundary conditions used, more delicate 
properties, like the droplet shape, are. Therefore in this section we re-derive the 
results of the paper [B] in the case of, more realistic, periodic boundary conditions. 

The calculation of free energy, expected values, and correlation functions for the 
spherical lattice gas is reduced, in a routine fashion, to calculation of the large-n 
asymptotics of an integral, see [6]. Introduction of new integration variables yj, 
j G V n in Eq. (JSJ) via the orthogonal transformation 



Xj 



J2wfy k , j G V n , (6) 

kev„ 



where the eigenvectors {wj,j G V n } are given by Eq. PJ, diagonalises the inter- 
action matrix. As a result, we obtain the following expression for the partition 
function 



9n(p) 



/oo roo I \ 

. . . / exp \/3J J2 hvl + P J2 VkVk fin(dy), 
J -°° \ kav n keVn j 



(k) 

where (p k = J2jev„ hjw) and 



]2 n (dy) = 6(y {hl _ l) VN-pN)5lY,y]- N ) II d Vr 

\jeVn j jev n 



Integration over 1/(1,1,. yields 

©n(p) = 



-7= / •••/ exp \pJdp 2 N + pjY, hyl + PpJ2 h i+PJ2 fkVk > 
ViV^-00 j-00 y k£Vn j&Vn k£Vn j 

8\T.'yi- m - p 2 )) n 'dy k , 

\kev n j keVn 

where primes indicate that summations /pro ducts do not include k = (1, 1, ... , 1). 
The integral representation for the delta function 

(, \ 1 r+ioc 

S„ vl - N(1 - p) )^^L dsexp 



s [N(i-fJ')-Y.'yl 



allows one to perform integration over the variables y k . However, we can switch the 
order of integration over the variables y k and s only after a shift of the integration 
contour for s to the right. The shift must assure that the real part of the quadratic 
form involving the variables y k , is negatively defined. On switching the integration 
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order, integrating over y^, and introducing a new integration variable z via s = ftJz 
one obtains 



Q n (p) = -^-= exp ftJdp 2 N + p p J2h j }[ — ) / dz exp [JV/3$ n (z)] , 

(7) 

where 

and c > d is the shift of the integration contour mentioned above. 

Depending on the value of ft, the large- n asymptotics of the integral ((71) can be 
found either by the saddle-point method, or by direct integration after an appro- 
priate change of variables. The saddle point of the integrand is a solution of the 
equation 

If the field {hj,j G Z d } is of the type ©, then 



1 v <f% < n-^b* 



Hence, for any z > d, as n — > oo, the sequence of the derivatives & n (z) converges to 

#{z) = J(l-p 2 )-±W d (z), 

where 

J—n J -it Z — 2^=1 COS0J v l=l ^7T 

is the Watson function. 

The function increases monotonically with z on [d, oo), and the location of 
its zeroes depends on the dimension d of the lattice. Namely, if d — 1,2, then the 
function $>'(z) has exactly one zero on the interval [d, oo) at a point z*(p) > d, for 
any ft > 0. If d > 3, then there exists a critical value 

of the parameter ft. If ft < ft c (p), then the function $'(V) still has exactly one zero 
on the interval [d, oo) at a point > d. While if ft > ft c {p), then the function $>'(z) 
is strictly positive on the interval [d, oo). 

In this section our aim is to investigate the natural state of the spherical lattice 
gas, which, as we shall see shortly, happens to be a mixed phase. Therefore we do 
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not utilize any symmetry-breaking perturbations and set hj = for all j G Z d untill 
Section 6. 

Application of the saddle-point method to the integral (J7]) is fairly straightfor- 
ward when there exists a saddle point z*(p) greater than d. In this case 

s*(z*W) 

and the saddle-point method yields 

-/„ = In 0„(p) = -L In ^ + J dp 2 + + O (n^lnn 

as n — >■ oo, where 

d ^ d dut 



L d (z) = J ...J Iniz-J^cosuAfl 



u=l J 1=1 27T 

When /3 > /3 C , the function $> n (z) still attains its minimum on the interval (A s , oo) 
at a point z* > A s , where A s = d — 1 + cos — is the second-largest eigenvalue of 
the interaction matrix T. However, the sequence of saddle points z* n approaches the 
pole of the integrand at z = A s , and the application of the saddle-point method 
becomes a bit more tricky. 

To find the asymptotics of the free energy we have to introduce a new integration 
variable ( via z = A s + (N^ 1 . The large-n asymptotics of the sum 

L { d N) (A B + CJV- 1 ) S 4E , ^(A. + CJV- 1 -A, 

is given by 

Lf (A, + CJV-) = H^fM + i £ " In (A, - X t ) - X E + o(iV - I ), 

where the double prime indicates that the summation does not include both the 
largest and the (2rf-times degenerate) second-largest eigenvalues. 
Hence 



n (p) ~ exp 



Npj(dp 2 + A s (l - p 2 )) -i]T"m (A s - A A 



x(iV-l)/2 J[(1 _ p 2 )/3 _ /3c]c 

rd-3/2 / \ 1 / e 



X 



"^{pj) mJo c 5 " c (10) 

where the integration contour C runs below the negative semi-axis (— oo, 0] — z 0, en- 
circles counterclockwise, and returns to — oo above the negative semi-axis (— oo, 0] + 
zO. 
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We see that the remaining integration does not contain the large parameter N 
any longer, and calculating the residue of the integrand at ( = we obtain 



1 

27ii Jc 



3 J[(l-p 2 )/3-/3 c ]C 



J' 



d-i 



(d-iy. 



(i - P 2 )/3 - t 



d-l 



We can use the same method to find the joint characteristic function s) of 
random variables Xj and xi. First, we express x(*> s ) = (exp(itXj + isxi)) n as the 
following integral 



PJ ( 71 \ 



(N-l)/2 



+ioo+e 



joo+c 



dz exp 



(11) 

Next, introduce a new integration variable ( via z = A s + 1 , and use the asymp- 
totic expansion for the partition function O n (p) to obtain 



X(t, s) ~ exp 
(d-1)! 



i(t + s)p - 



t 2 + s 2 
4/3JN 



E 



A« — At 



E 

fcev„ 



(fe) (fc) 

A a — At 



X 



2m 



7c exp [ c " I 1 " ^ J V ( s } S cos — ~ 



d( 

(12) 



The remaining integral can be expressed in terms of the Bessel function 



2m' Jc 



exp 



c- 



d( 



c 



m+1 ' 



however, calculation of the moments can be done easily by differentiation under the 
integral sign in Eq. ( 112"]) . 

Setting s = and passing to the limit iV — > oo we see that the random variables 
Xj of the spherical model on the infinite lattice have the common characteristic 
function 



X(t) = exp 



itp — 



t 2 



4/3 J 



W d {d) 



(d-l)\ 
2tt% Jc 



exp 



c 



1-p 2 

2C 



1 - 



c 1 .2 



(13) 



Differentiation over i shows that the expected values are given by (xj) = p, and the 
variances are given by 



Var(a;j 



W d (d) 1-p 2 f _ & 



2/3J 



11 



More importantly, differentiation of the joint characteristic function shows that 



Covfo, x t ) = (l-p 2 ) (l-^Y-T, cos ^ Q 



(3 d 



u=l 



n 



Note that Cov(xj,xi) -fa 0, as dist (j, X) — > oo, for any /3 > (3 C . Therefore, the 
lattice gas is in a mixed (or degenerate) state, and in order to calculate the values 
of observable quantities we have to single out pure phases. 

In the case of the stretched lattice ([1]) the second-largest eigenvalue of the inter- 
action matrix is only twice degenerate, and we obtain 



X(t, s) 



exp 



i(t + s)p - 



t 2 + s 2 
A/3JN 



fcex n 



A s — A j 



ts 

2/3J 



E 



(fc) (fc) 
A s — At. 



Jo 



\ 



(14) 



where Jq{z) is the Bessel function of order zero. Hence, for 1 <C \j\ — h\ <C n we 
have 



X{t,s) ~ £exp 



i(t + s) 



\ 



2(1 -p 2 ) fl-^J sin27rt/ 



P. 

w d (dY 

2/3J 



where the random variable U is uniformly distributed on the interval [0, 1], and it 
is independent of the normal random variables TV} and A/j. 

The last expression admits the following probabilistic interpretation. Let 



j/n 



3l 32 

n n 



3d 
n 



) ->■ 7 = (7i,72,-- • ,7d) 



7 d J, as n ^ oo, 



then in the continuous limit the field of random variables Xj converges to the random 
function 

W d (d)\ 



\ 



2(l-P 2 )(l-|)sin2 I rU 7 +A^U ^ 



where 7V 7 (a, b) is a continuous function of independent normal random variables 
defined on the rectangle 

< 71 < 1 + < 72 < 1, . . . , < 7d < 1. 

The random variables £7 7 are uniformly distributed on the interval [0,1], any pair 
t/ 7 , U a becomes perfectly correlated as 71 —> a\. While if 71 7^ a±, the random 
variables U 1 , U a are neither independent, nor perfectly correlated. 
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5 Droplet shape at zero temperature. 



As it often happens, the situation at zero temperature is simpler than in the case 
> 0. If = 0, then the pure states of the spherical gas are configurations 
x 7v = {xj, j G V n } which minimize the energy H n subject to the spherical and the 
density constraints. The corresponding minimization problem is solved using the 
Lagrange function 

C({xj,j G V n }, a, b) = - J £ T ikXjX k + a £ x) - N J + bl £ Xj - pN 

where a and b are Lagrange multipliers. Introducing the variables {yj,j G V n }, see 
Eq. (jfJJ), we get rid of the density constraint, and write down the Lagrange function 
in the following diagonal form 

£({ yj ,j G V n }, a) = - JJV^Ad,!,...,!) - J £ \ 3 y) + a ( £ ' y] - N(l - p 2 ) ] , 

where the primes indicate that the summations do not involve j — (1,1, ... , 1). 

Differentiating over {yj, j G V n \ (1, 1, . . . , 1)} we obtain the following system of 
equations for stationary points 

2J\ jVj = 2a Vj , j G V n \ (1, 1, . . . , 1). 

If a 7^ JAfc, then the only solution of the system, yj = for j G V n \ (1, 1, . . . , 1), 
violates the spherical constraint. Therefore to solve the constraint minimization 
problem we have to set the value of the Lagrange multiplier to a = JXk, for some 
k G V n \ (1, 1, . . . , 1). Using the spherical constraint we obtain for such a value of a, 
that 

]T y] = N(l-p 2 ), and Vj = 0, if JX^a. 

j:JXj=a 

The global minimum of — X^ e y„ ^jVj ^ s obtained when a = JA S , where A s is the 
2c? times degenerate second largest eigenvalue of the interaction matrix. Therefore 
the configuration {yj,j G V n } minimizing — Y^j^Vn^jV] subject to the spherical and 
density constraints has the following components. The component yt l ^ = p\fN , 
due to the density constraint. The 2d components y*, with j such that Xj — X s , 
are the coordinates an arbitrary point on a 2 d- dimensional sphere with the center 



at the origin and the radius J N(l — p 2 ). All the remaining N — 2d — 1 components 
of the configuration {y*j,j G V n } are zeroes. 

Going back to the variables {xj,j G V n } we obtain 



Xj = w 



(1 ' W) pv/iV+ £ wfyl = p+ £ wfyt. 

k:\ k =\ s k:\ k =\ a 
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The 2d orthonormal eigenvectors corresponding to the eigenvalue A s are given by 

~2 



w 



± 



cos 



27X{J V - 1) ± 7T 

4 



n 



1,2,. ..,d. 



Let us denote yj^ the component yl such that = w^. Then the components of 
pure zero-temperature states x^ are given by 



x 



j = p + V^E cos 

d 



~~ 4 



n 



= p+y^EV(^) 2 + (^) 2 co6 

where the phase shifts a u are such that 



V v + cos 
27r(> - 1) 



n 4 



(15) 



cos a v 



IT 



v / (^) 2 + (^ + ) 2 ' 



and sin a 



v / (^) 2 + (^ + ) 2 ' 



The phase shifts {a u }f, =1 merely set the "center" of the pure state Xn, which 
we are free to choose at will because of the translation invariance of the Hamil- 



tonian of the spherical lattice gas. The multipliers y (y„) 2 + iyt) 2 play a much 
more important role, they determine asymmetry of the droplet across dimensions. 
If (y^) 2 + {yt) 2 — N(l — p 2 ), and yf = 0, for v ^ k, then the states are maximally 
asymmetric. Namely, they have the cosine shape in the dimension v — k, and they 
are translation invariant across other dimensions, see Fig. 1, 



p + v /2(l-p 2 )cos 



Mjk - 1) 



n 



The opposite extreme is the symmetric case (y v ) 2 + (y„) 2 = N(l — p 2 )/d, for 
v — 1,2, ... ,d, when the droplet has a rounded-square shape, see Fig. 1, 



Xj — p + 



'2(1 -p 2 )* 

Ecos 



d 



2n(j v - 1) 



a. 



n 



(16) 



In intermediate situations, the droplet has a kind of elliptic shape, see Fig. 2. 

It turns out that the above zero-temperature pure (macro) states are stable with 
respect to heating up of the system to a temperature not exceeding the critical value 
Pure (macro) states for > are in one-to-one correspondence with zero- 
temperature pure states. In the next section we single out pure states for [3~ l > 
by quasi-averaging. 
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Figure 1: The shape of the maximally asymmetric and symmetric droplets. 




Figure 2: A droplet of an intermediate shape. 

6 Quasi- averages. 

Pure (micro) states of thermodynamic systems can be obtained with a help of quasi- 
averaging pp. The basic recipe of quasi- averaging looks like this: switch on an 
external field (of magnitude e) removing the degeneracy of equilibrium state, pass 
to the thermodynamic limit, switch off the field by sending its magnitude e to zero 
over positive values {e \. 0). It turns out that in the case of the spherical lattice gas 
this procedure often deforms macro states (droplet shape). Therefore, we have to 
use a more potent version of quasi-averaging [2J [3], where the field magnitude is a 
decreasing function of the volume N. That is, the symmetry-breaking external field 
is switched off together with the thermodynamic limit. 

Motivated by the results of the previous section, we begin with the Hamiltonian 
fl2l) where the external field 



:V2 



■ cos 



v. 



« 



?? 



3 = Ui,32,---Jd) e v n 



(17) 
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is proportional to an eigenvector corresponding to the second-largest eigenvalue, A s , 
of the interaction matrix. The dimensionality of the eigenspace corresponding to A s 
is 2d, and we are free to choose the value of the parameter a G [0, 2ir), and instead 
of ji we can use j u , with v = 2, 3, . . . , d. Different choices of a and v yield different 
pure macro states. 

For the characteristic function of random variables Xj and xi one obtains an 
integral similar to Eq. ( 11 ip . 



X{t,s) 



exp(/3Jdp 2 N)e itp+isp 0J ( n \ 
~2m \Jj) 



(N-l)/2 



&n(p)VN 



X 



+ioo+c 



— ioo+c 



dz exp [Nj3Q> n (z, e)] exp 



keVr 



W(z - A fc ) 



(18) 



where 

Np$Jz,e) = NP$Jz) + 



1 £ 



sw 



(2,1,..,1) 
I 



2J n s-d/2 



z — A t 



+ 



fls 2 



AJn 25 - d {z- A s ) 

=1pn(z) 



The innocent-looking term ip n (z) makes a huge difference for the evaluation and 
behaviour of the integral ([181) when (3 > f3 c . Indeed, in the scale z = A s + (n~ 5 
the magnitudes of variation of N(3Q n (z) and ip n {z) with ( become comparable. 
Moreover, if S < d, then ip n (z) prevents the saddle point 2* from approaching the 
pole at A s any closer than the distance 0{n~ 5 ). 

Introducing a new integration variable ( via z = A s + (n~ s one obtains 



N^ n (X s + (n~ 5 , e) = iV/3$„(A s + (n 



-S\ 



eV2N \ 
i— — — { t cos 



2JC 



+ s cos 



2tt(/i - 1) 



a 



n 



n d ~ d [3e' 



Gathering the terms of the order n d 5 and differentiating over £ we obtain the 
following equation for the saddle-point (*: 



(l-p 2 )J(/3-/3 c ) 



(3e 2 
4JC 2 



0. 



Hence the saddle point of $> n (z, e) is located at 
* _ \ |e| 1 

Z n — A e + 



2Jn l 



y J(l-p 2 )(l-(3 c /(3) 



+ o(n 



Application of the saddle-point method to the integral ( fl8|) and to the analo- 
gous integral for the partition function Q n {p) shows that the large-n asymptotics of 
x{t, s) coincides with the characteristic function of a bivariate normal distribution 
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B(fii, Hi\ Vi, V2, c). That is, for large n the joint distribution of random variables Xj 
and Xi is asymptotically normal with mean values /ii and /i2, variances v\ and u|, 
respectively, and covariance c. Under the assumption that the first components of 
the nodes j and I scale with n as ji ~ •yn and Zi ~ An, the expected values of Xj 
and xi are given by 



Mi = p + v /2(l-p 2 )(l-/3 c //3) cos (2vr 7 - a) , (19) 

/i 2 = p + y/2(l -p 2 )(l -/3 c //3) cos (2ttA - a) . (20) 

Note that, provided 5 G (0, d), the expected values (droplet shape) do not depend 
on the magnitude of the external field en~ s . Therefore it is reasonable to conclude 
that the symmetry-breaking field (fl7|) selects one of the possible droplet shapes, but 
it does not deform the droplet. 

The variances vf, v\ and covariance c(j,l) are given by the usual formulae for 
the pure phases of the spherical model below the critical temperature 

2 2 W d (d) 
Vl = V 9 = , 

1 2 2/3J ' 

exp 



__ , , -: [i?Z = iUu - lu)u>v] «cU T(d/2 - 1) 
C[J ' L) 2/3J J-,--]-, d-Zi =l cosu v 2tt ~4/3Jtt^7 2 ' 1 J 

where r^y is the Euclidean distance between j and Z. 

The equations (fl9~j) and fTSOj) describe the shape of a localised droplet obtained 
with a help of the (generalized, see [21 E]) method of quasi-averages. To double 
check that the droplet shape (fT9l is not deformed by the external field we would 
like to examine the sensitivity of the droplet shape to the type of external field used 
for quasi-averaging. For this purpose we now repeat the above calculations for a 
technically more demanding case 

h=i £n ~ S ' if -? 1 = 1 ' n e v (22) 

h i \ 0, otherwise, 3 ^ Vn ' [ZZ) 

where e > 0, and 5 > 0. First, we have to find the range of values for S, such that the 
field (1221) is strong enough to fix the location of the droplet of condensed "spherical 
matter", but, at the same time, it is weak enough not to deform the droplet shape. 

The partition function Q n (p) is still given by Eq. (j7j), but the coefficients ip^, 
k eV n are now given by 

^ = £ n^ w n^,l). 

In order to find the location of the saddle-point of the integrand in Eq. (J7]) we have 
to first investigate the behaviour of the sum 

<z( z \ = J_ v ; Vk = v - 

V ; 4JiV fc ^ n z - A fe 4Jn 2 ( 1+<5 ) ^ 2 z - d + 1 - cos [2tt(A; 1 - l)/n] 
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in the vicinity of the point z = d. Fortunately, this sum can be calculated exactly, 
see [ID], 



4J n 2(l+S) 



2n 



xMz) + 1 



%2{z) — x\(z) x'^z) — 1 z — d 

Hence, for any fixed z > d, the (derivative of the) sum S(z) = 0(n~ 1 ~ 2<5 ) produces 
only a vanishing contribution to the saddle-point equation (jSJ). 

The contribution of the sum £(z) becomes significant below the critical temper- 
ature, where we have to introduce a new integration variable ( via z = A s + (n~' y 
before application of the saddle-point method. In order to find the right rescaling 
for the integration variable z (which, as we shall see shortly, depends on the value 
of 5 in Eq. (1221) ). we have to analyse the behavior of £(z) for different values of 7. 

If < 7 < 2, then 



£(A S + Cn- 7 ) ~ 

If 7 = 2, then 
S(A S + C^ 7 ) 



n 



-1-2(5+7/2 



as n — y 00. 



e 2 n- 25 



COth J(/2 - 7T 2 



U 



C/2 



7T Z 



C/2 - n 2 



as n — > 00. 



Note that the r.h.s. of the last equation does not actually have a singularity at 
C = 2tt 2 . Therefore to find £(A S + C n ~ 7 ) f° r C < 27r 2 one can use the analytic 
continuation 



S(A S + Cn" 



£ 2 n -28 I cot ^7T 2 - C/2 



8J 



7T 2 - C/2 C/2 



as n — y 00. 



Ti 



Finally, if 7 > 2, then the main asymptotics of S(A S + C^ 7 ) comes entirely from 
the first and the last terms of this sum, 



S(A S + C^ 7 ) 



£ 2 n7 -2(l+5) 



2JC 



as n — y 00. 



Under the same rescaling z = A s + (n 7 the first two terms of <& n (z) become 
J(A S + Cn^)(l - p 2 ) - £ ' H^s + C^ 7 - A*) ~ 



JA s (l-p 2 



L d {d) 

2/3 



+ r^ 7 C 



J(l-p 2 )- 



2/3 



+ n' 



- 2 7r 2 iy,(d) 

/3 ' 



Hence, if < 5 < 1, then the contributions of all C-dependent terms in $ n (A s + C^ 7 ) 
are of the same order when 7 = 2(1 + 28)/3, and the saddle-point equation for C is 
given by 

2 w„ n , (3 e 2 



J(l-p 2 )(/3-/3 c ) 



U2V2C 3 / 2 ' 
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The positive solution of this equation is given by 



\4/3 



2 \2J) [(1 -p*)]*/ 8 " 



(23) 



If 8 = 1, then we have to set 7 = 2, which yields the following saddle-point 
equation for £ 



J(l-p 2 )(/3-/3 c ) 



(3 £ 2 d /coth JC/2 



7T 



1 



C/2 - 



(24) 



Finally, ifl<5<d — 1, then we have to set 7 = 1 + 5, and the saddle-point 
equation for £ is given by 



(3e 2 
2JC 



The positive solution of the above equation is given by 

e 



Jj2(l-f3 c //3)(l-p*) 



Now that we know the behavior of the saddle point z* n of the integrand in Eq. ([7]) 
we can find the thermodynamic limits of various macro- and microscopic quantities. 
But before we are able to apply the saddle-point method and find the characteristic 
function of an arbitrary pair of random variables Xj and x\ we have to calculate the 
sum 



E 



(k) 



z — Xi 



appearing in Eq. ( ITT]) . After some elementary transformations the above sum reduces 
to 



n 



e_ " cos [27r(ji - l)(fci - l)/n] 

1+5 2^ 



:2 z — d + 1 — cos [27r(fci — l)/n] 
and the summation technique from [10] yields 

2s 



I r i- ii+i 

2# - 1 



n 1+5 (z-d)' 



On application of the saddle-point method to the integral in Eq. ( II ip one obtains 
the following expression for the characteristic functions of random variables Xj and 

Xf. 



X{t,s) 



exp 



4/3J(z*-\ k ) 



+ it (E, «)+p)+i S (E,«) + 



• (25) 



19 



Hence, as it is usually the case for pure states of the spherical model, for large values 
of n the random variables {xj,j G V n } have asymptotically normal distributions. 

If 5 G (0, 1), then the saddle point is given by z* n = \ s + (*n~ 2(1+2<5)/3 , see Eq. 
( 12"5|) , and the large- n asymptotics of the expected values of the microscopic variables 
(the multipliers of it and is in Eq. fl25]) ) are given by 

JlS)/3 r 



( X j) 



en y 



h2<5)/3 



-A 



2JV2C 

if ji < n/2. In this case the droplet shape is clearly deformed by the external 
field, since virtually all "spherical matter" gathers in a narrow layer of the width 
~ ra ( 1 + 2,5 )/3 ar ound the hyperplane j\ = 1, where the field is applied. 

If 5 = 1, then the saddle point is given by z* n = X s + (*n~ 2 , where is a solution 
of Eq. ([21]), and the large-n asymptotics of the expected values of the microscopic 
variables are given by 



cosh 



(l-2 7l )^C72-7T s 



1 



2(C* - 2tt 2 ) sinh J(*/2 - vr 2 C 



2tt 2 



(26) 



if j l ~ ry lTl _ j n this case i}yQ droplet shape is still deformed by the external field, 
since it depends on the field's magnitude e. 

Finally, if 5 G (1, d — 1), then the large-n asymptotics of the expected values of 
the thermodynamic variables are given by 



( Xj ) ~ p + V / 2(l-p 2 )(l-/3 c //3) cos (2vr 7l ) . 

if j 1 ry-^n. In this case, apparently, the droplet shape is not deformed by the 
external field, since it does not depend on the field's magnitude e, and it is the same 
as in Eqs. (fT9|) and (|20|). 

Thus, the droplet shape is sensitive to the magnitude of the symmetry-breaking 
field only if it scales with n as n~ s with 5 G [0,1]. If the field is switched off 
reasonably fast, 5 G (1, d — 1) in Eq. ( 1221) . the droplet shape is not deformed by the 
field (does not depend on the magnitude of the field). If the field is switched off 
too fast, S > d — 1, then it does not transfer the thermodynamic system into a pure 
state. 

It is appropriate to stress at this point that different configurations of the sym- 
metry breaking field {hj : j G V n } can produce different pure macroscopic phases. 
Therefore the droplet shape is sensitive to the type of field used for quasi-averaging. 
The eigenvector fflTl) is only one of 2d orthogonal vectors spanning the eigenspace 
corresponding to the second-largest eigenvalue of the interaction matrix. In this 
respect, the symmetry breaking field ( J22|) is special, because it is orthogonal to all 
of the remaining 2d — 1 eigenvectors corresponding to the eigenvalue X s . If instead 
of the field ff22l we decide to use a configuration {hj : j G V n }, which has a non-zero 
scalar product, say, with 



cos 



Mh - 1) 



— a 



n 



j = (ji,32,---,jd) e V n 
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then the corresponding vector of expected values (xj), j e V„ will contain a compo- 
nent proportional to cos(27rj 2 /n), cf. Eq. ([16]) . 

7 Discussion. 

The investigation of droplet shape performed in this paper demonstrates the impor- 
tance of decomposing mixed states into constituent pure components. It appears 
that, often, mixed states are mathematical constructions reflecting absence of suf- 
ficient information about the system under investigation. A priori we do not know 
which of many possible pure states of the system will be actually observed in an ex- 
periment. As a consequence, a direct mathematical solution of the model produces 
a mixed state. At the same time, if an experimentalist actually performs measure- 
ment in the corresponding physical system, the obtained results would be described 
by one of pure states. 

Whether this point of view is valid or not, depends, of course, on how quickly the 
system under investigation can transit from one pure state into another. Namely, if 
the typical transition time is much larger than the typical experimental observation 
time, then measurements yield results corresponding to one of pure states. If the 
observation time is much longer than the transition time, then the observed values 
are likely to correspond to a mixed state. 

Investigation of dynamical properties of the mean spherical model were per- 
formed recently, see [5]. It is possible to calculate the interstate transition times 
using the technique proposed there. We have not attempted yet those calculations, 
and we do not even know if the droplet shape in the mean spherical and the con- 
ventional spherical model are similar. But, it does not seem inconceivable, that the 
transition time from one pure state into another is much greater than the typical 
measurement time. Indeed, investigations of large deviation probabilities show that 
in order to transit from one pure state into another within the conventional spherical 
model one has to overcome a free-energy barrier that grows exponentially with n d ~ 2 . 

It is well known that gaussian approximations correspond to leading orders of 
low-temperature expansions of various continuous models. Therefore, one can hope 
that the behaviour of realistic continuous models, for instance 0(n) models, is quali- 
tatively similar to the findings of the present paper. On the other hand, the common 
knowledge is that physical substances in liquid phases form droplets which bound- 
aries are sharp in the macroscopic scale. Most likely, the sharpness of boundaries 
should be attributed to the atomistic (discrete) microscopic structure of physical 
substances. Therefore, in this respect, discrete lattice gas models of Ising type 
provide much more realistic description of various substances in the liquid state. 
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